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In this Letter, we propose a new generalized Ricci dark energy (NGR) model to unify Ricci dark 
energy (RDE) and XCDM. Our model can distinguish between RDE and XCDM by introducing 
a parameter /3 called weight factor. When 13 — 1, NGR model becomes the usual RDE model. 
^SJ ' The XCDM model is corresponding to /3 = 0. Moreover, NGR model permits the situation where 

neither /3 = 1 nor j3 — 0. We then perform a statefinder analysis on NGR model to see how /3 
04 ' effects the trajectory on the r — s plane. In order to know the value of j3, we constrain NGR 

model with latest observations including type la supernovae (SNe la) from Union2 set (557 data), 

D ■ 

' baryonic acoustic oscillation (BAO) observation from the spectroscopic Sloan Digital Sky Survey 

(SDSS) data release 7 (DR7) galaxy sample and cosmic microwave background (CMB) observation 
CN ■ from the 7 -year Wilkinson Microwave Anisotropy Probe (WMAP7) results. With Markov Chain 

Monte Carlo (MCMC) method, the constraint result is l3=0.08to.li{^o-)to.2U'^^)^ which manifests 
the observations prefer a XCDM universe rather than RDE model. It seems RDE model is ruled out 
in NGR scenario within 2a regions. Furthermore, we compare it with some of successful cosmological 
models using AIC information criterion. NGR model seems to be a good choice for describing the 
universe. 
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1. Introduction 

^: . . . ... . n . . 

, Various cosmic observations suggest our universe is undergoing an accelerated expansion 1". To explain this 
phenomenon, people introduce an exotic component with negative pressure known as dark energy. The simplest 
dark energy model is cosmological constant (A) 2] or XCDM model where dark energy has an arbitrary equation 
. 5^ ' of state (EOS) lux- It fits all kinds of observational data well while it is confronted with theoretical problems such 
>< , 

. . .' proposed including quintessence quintom [5], phantom fl, GCG R and so on. In principle, dark energy is related 
to quantum gravity Q]. But until now, a self-consistent quantum gravity theory has not established. Nevertheless, 
the holographic principle Q is thought to be a reflection of quantum gravity. Motivated by this, holographic dark 
energy has been proposed. It embodies the relation between UV cut-off and IR cut off. However, how to choose the 



(N 



IR cut-off is a problem. Cohen et al. pLfl] first chose Hubble scale as IR cut-off. Hsu and Li ll| pointed out it can not 
give an acceleration solution. Li then suggested the future event horizon_as IR cut-off 12|. Basing on causality, Cai 



proposed agegraphic dark energy [13| and new agegraphic dark energy jlil. Furthermore, Gao et al. [3l proposed a 
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holographic dark energy from Ricci scalar curvature. In RDE model, the IR cut-off is determined by a local quantity. 

Nowadays, all the models above seem to be consistent with current observations. Usually, we estimate models 
through the or information criteria like BIC and AIC [igI. In this Letter, we find XCDM model and RDE model 
can be related by a parameter /3, thus we can estimate them through constraining l3. The distribution of /? can reflect 
which model is better. For example, if the best-fit value of /3 is close to 1 and is not within 2-a range, we can say 
the observations support RDE model rather than XCDM model. We now give some similar examples. In order to 
know whether ACDM is right, people free the EOS parameter and constrain it with observations. If the result is close 
to -1, we can say ACDM is still a good choice. However, if the EOS parameter tends to -2, then ACDM should be 
suspected. Likewise, for purely dimensional reasons, Granda and Oliveros 17[ proposed a new IR cut-off. Wang and 
Xu [3| give the constraint results which suggest the coefficient of is two times larger than the one of H, thus 
ruling out the SRDE model 



19| . In RDE model, the density of dark energy is proportional to Ricci scalar or the sum 
of traces of energy-momentum tensors of each component. Since the trace of radiation is 0, we can ignore its impact 
on space-time curvature. RDE model suggests the weights of dark energy and matter are the same, while XCDM 
model suggests only the trace of dark energy can affect its density. Therefore, what on earth is the weight of matter 
(0.5 or 0?) is an interesting thing we want to know. Motivated by this, we free the weight of matter as an arbitrary 
parameter called weight factor. 

The rest of the Letter is organized as follows. In Section 2, we give the dynamics of the new generalized Ricci dark 
energy model. In Section 3, we give a statefinder diagnostic. In Section 4, we introduce the observational data we use. 
The constraint results are shown in Section 5. At last, we give the discussion and conclusion in section 6. Throughout 
the Letter, the unit with light velocity c — 1 is used. 



2. New generalized Ricci dark energy model 

We assume the universe is flat and described by Friedmann-Robertson- Walker (FRW) metric. For RDE model, the 
density of dark energy is proportional to Ricci scalar 

Pde oc R. (1) 
Considering Einstein field equation can be expressed as 

R = SttGT, (2) 
where G is Newtonian constant, T is the sum of traces of each component, RDE model can be expressed as 

Pde oc Tde + Tm OC Pde - + Pm- (3) 

From this equation, the coefficients of Tde and Tm are both 1, which means the weights of dark energy and matter 
are the same. We now change the weight of matter, the equation becomes 

Pde = a{Tde + PTyn) = a(Pde ^ ip + P Pm) , (4) 

/3 here is the weight factor we introduce which reflects the relative weight of matter to dark energy. If /? = 1, it 
becomes the usual RDE model. When /3 = 

Pde = a{pde - ip), (5) 
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equivalently, 

Pde oc p, (6) 
the NGR model becomes XCDM model. For simplicity, we define dimensionless quantities = — , ^de = — , 

Pc Pc 

where pc = |^ is the critical density of the universe. H is Hubble parameter, subscript "0" represents the quantity 
today. The Friedmann equation can be expressed as 

E'^ — ^de + ^m, (7) 

where E = -jj-- 

-Ho 

The energy-momentum conservation equation can be expressed as 

+ 3i/(l + Wi)!^, = 0, (8) 

subscript "i" represents dark energy or matter. Then we get 

^, _ (4a- l)J7rfe + a/3a„o(l + 2)^ 



where fi'^^ = With the initial condition 



a{l + z) 



fideO + ^mO = 1, (10) 



we can obtain the evolution of Vide with respect to redshift z 



1 — a 1 ~ a 



The EOS parameter can be obtained by 

lij^^ = — 1 -A- (^ ^ •/' 

i^d. 

and the deceleration parameter 



a;de = -l + (1 + ^)5^, (12) 



.= ^1 + 1^). (13) 

In order to exhibit the effects of /3, we fix J7,„o = 0.27 and i>Jdca = — 1 and plot the evolutions of lOdeiz)^ Hubble 
parameter H{z) and density parameters defined as Vti/E"^ in Fig. [T]and Fig. [51 

3. Stateflnder diagnostic 

Statefinder diagnostic is a useful method to differentiate effective cosmological models since these models are all 



seen to be consist with current observations. It was first introduced by Sahni et al. [20[. This method probes the 
expansion dynamics of the universe through high derivatives of scale factor a . The dimensionless statefinder pair 
{r, s} is defined as 

_ 'a _ r — 1 
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The NGR Model 


Q n 




and v"^ - 


13 = 0.5 

13 = 1 


0.284ini^(la)lHi^(2a) 
0.280lHf9(l-)lHi?(2-) 
0.287lH^3t(l-)lHf4(2a) 
0.296lH^3^(l-)lHrr(2-) 


0.235lHf9(l-)lHf3(2-) 
0.2461H?6(1-)1H3?(2-) 
0.195lHf9(l-)lH2?(2-) 
0.16llHf4(l-)tHS(2-) 


0.081^2? (l'^)tS;i(2'^) 531.710 

* 532.238 

* 539.734 

* 558.834 





FIG. 2: The evolutions of the Hubble parameter in units of Hacdm{z) (left) and the density parameters (right), flmo = 0.27, 

I^deO = — 1. 



Since the scale factor depends on the space-time manifold, the statefinder is a geometrical diagnostic. Different 
models are corresponding to different trajectories on the r — s plane. For example, the spatially flat ACDM model 
are corresponding to a fixed point on the plane, 



{s,r} 



= {0,1} 



(15) 



ACDM 

Statefinder has been applied to various dark energy models including quintessence, quintom, GCG, braneworld 
model and so on 2l|. We now turn to statefindr diagnostic for NGR model and find the effects of /3. The statefinder 
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parameters can also be expressed in terms of the total energy density and the total pressure 



r = 1 



^>iptot+p) P 



"^Ptot Ptot 

where we ignore the pressure of matter. 

Combined with the dynamics we discussed in section 2, we have 



{ptot +p) P 

P Ptot 



(16) 



r!,„o(l + + ^i^^ 3/3f]„,o(l + zf + (1 - i)n'de 



and 



, 9 -nde 

r = 1 H 



3 



(17) 



(18) 



2 + ^mo{l + Z)^ 

In order to plot the statefinder plane, we fix the current EOS of dark energy and the density of matter as tudeo — ~^ 
and flmo ~ 0.27, respectively. 

In Fig. [3l we can see with the increase of the value of /3, the corresponding s becomes smaller, and the range of 
the trajectory becomes larger. For XCDM model, we choose the initial condition as uJdeQ — ~1, it is regarded as 
ACDM model here. The dots represent the points today which are linear to fl. r=0.865, 1.135, 1.27 and 1.405 for 
/3 = —1/3, 1/3, 2/3, 1, respectively. Our results are consist with the RDE case 
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FIG. 3: The r - s plane for NGR model with = -1/3, 0, 1/3, 2/3, 1, respectively. 
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Current observational data 



4.1. Type la supernovae 



SNe la has been an important tool for probing the nature of the universe since it first revealed the acceleration of 
the universe. The current data (Union2) is given by the Supernova Cosmology Project (SCP) collaboration including 
557 samples [2^. The distance modules can be expressed as 

/i = 51og(dL/Mpc) + 25 , (19) 

where cIl is the luminosity distance. In a flat universe, it is related to redshift which is a observational quantity 

= (1 + z) / dz'/H[z'). (20) 



"'0 

We choose the marginalized nuisance parameter |2_4| for y^: 

xInc = A-^, (21) 

where A = X),^^ (m'^''*'' - M*'"'°''^)Vcr2, B = J^T'^ (m'^''*'' " m'^™'^)/^! , C = ^ l/a^ a, is the Ict uncertainty of 
the observational data. 



4.2. Baryon acoustic oscillation 



For BAO, the distance scale is expressed as 25 1 



, X 1 r ^ / /"^ N2nl/3 , ^ 

h^IemI/.ew" - '''' 

and baryons were released from photons at the drag epoch. The corresponding redshift is give by 

I — I 

where &i = 0.313(^7mo/^^)""■'*^^[l + 0.607(17,„o/i^)°-^^1"^ and 62 = 0.238(r2,„o^^)°-^^^ [26]. For BAO observation, we 
choose the measurements of the distance radio (dz) at z = 0.2 and z = 0.35 t27|]. It can be defined as 



-'-B.,,r '''' 

where Tg{zij) is the comoving sound horizon. The SDSS data release 7 (DR7) galaxy sample gives the best- fit values 
of the data set (do.2, do. 35) [27 



p 

^ matrix — 


' d"o.2 






I do. 35 y 





0.1905 ±0.0061 
0.1097 ±0.0036 



(25) 



27| 



The value of this BAO observation from SDSS DR7 can be calculated as 

XbAO = ^PmatrixC'matrix ^APmatrix, (26) 

where APmatrix = Pmatrix — Pmatrix, and the Corresponding inverse covariance matrix is 

_^ I 30124 -17227 , 

Cinatrix =1 ■ (27) 

-17227 86977 
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4.3. Cosmic microwave background 



For CMB, the acoustic scale is related to the distance ratio. It can be expressed as 



rs(z*) 



(28) 



where ^^(z,) Hq Cg{z) / E{z)dz is the comoving sound horizon at photo-decoupling epoch. The redshift 

corresponding to the decoupling epoch of photons is given by 28 1 



= 1048[1 + 0.00124(r2b/i2)-o '^38(^ _,_ 5^(l7,„o^^)^')] 



where gi = 0.0783(i]6/i^)-"-^'^«(l + 39.5(f7h/i^)-"-7&^)-i, = 0.560(1 + 21.1(f7fc/i^) 
R is expressed as 



,2M.81N-1 



29| 



(29) 



The CMB shift parameter 



1/2 



dz 



(30) 



For the CMB data, we choose the data set including the the acoustic scale (la), the shift parameter (R), and the 
redshift of recombination (z»). The WMAP7 measurement gives the best-fit values of the data set j3o| 



CMB 



R 



/ 302.09 ±0.76 \ 

1.725 ±0.018 
V 1091.3 ±0.91 / 



(31) 



The value of the CMB observation can be calculated as 

XcMB = APcmbCcmb~^APcmb, 
where APcmb = Pcmb — Pcmb, and the corresponding inverse covariance matrix is 

/ 2.305 29.698 -1.333 \ 



(32) 



■^CMB 



29.698 6825.270 -113.180 
\ -1.333 -113.180 3.414 / 



(33) 



5. Constraint results 

We choose the common cosmic observations including SNe la, BAO and CMB to constrain the NGR model. We 
use the usual maximum likelihood method of fitting with Markov Chain Monte Carlo (MCMC) method. The code 
is based on CosmoMCMC 3l|- The total can be expressed as 

= xInc + xIao + XcMB • (34) 



We show the 1-D probability of each parameter {flmOi and /3) and 2-D plots for parameters between each other for 

035(1c^)-0.048(2ct), a = 0.235_o.o3g(lcr)_o.o53( 



the NGR model in Fig. H The constraint results arc n„M = 0.284t°:°i^^(lcr)t°;°^g(2cr), a = 0.235lJ5;^3^(lcr)l^;j;^^(2cr), 



/3=0.08^Q'2i(lo')to'28(2f). We can see /3 = is within l-cr range and /? = 1 is ruled out within 2a regions. Moreover, 
we further fix the value of /3 in three cases: /3 = (XCDM), /3 = 0.5 (the situation NGR model permits) and /3 = 1 
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FIG. 4: The 2-D regions and 1-D marginalized distribution with the 1-a and 2-a contours of parameters fimo, ce and /? in NGR 
model, for the data sets SNe+CMB+BAO. 

(RDE). The results are plotted in Fig. O Fig. Inland Fig. [3 respectively. Numerical results are shown in Table. H] 
We can see that when (3 becomes larger, the corresponding xLin becomes larger quickly. The xLin of R-DE model is 
558.834 while Xmin of XCDM model is only 532.238. We can also see when /3 becomes larger, the density of matter 
becomes larger and parameter a becomes smaller. Our constraint results are consistent with 



32|. 



6. Discussion and conclusion 

In this Letter, we propose a new generalized Ricci dark energy model based on the weight of matter. This model 
contains both Ricci dark energy model and XCDM model through weight factor (3. /3 = and (3 = 1 are corresponding 
to XCDM model and RDE model, repectively. Moreover, NGR model permits an arbitrary value of (3. If we fix the 
EOS parameter today uj^eo = —1 and ilmo = 0.27, which seems reasonable for all kind of observations, the larger (3 
becomes, the faster EOS parameter w^e tends to 0. Besides, deceleration parameter becomes smaller in the future, 
Hubble parameter becomes larger and density parameter of dark energy becomes larger. The observations can give 
us the distribution of [3, which provides a criterion for testing XCDM and RDE. It is similar to testing the distance- 
duality relation [3^. Both of them set the key parameter free. We use the latest observational data including SNe la, 
BAO and CMB to constrain NGR model. The constraint results tend to supporting XCDM model or even ACDM 
model (corresponding to /3 = and a = 0.25) rather than RDE model. We can conclude that RDE model is ruled out 
by the observations we select in NGR scenario within 2a regions. For future study on this problem, we hope more data 
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FIG. 5: The 2-D regions and 1-D marginalized distribution with the l-cr and 2-(t contours of parameters flmo and a in XCDM 
model, for the data sets SNe+CMB+BAO. 

Model Number of parameters Xmin AAIC 



ACDM 1 532.313 

XCDM 2 532.238 1.925 

RDE 2 558.834 28.521 

GCG 2 532.159 1.846 

GPL 3 531.804 3.491 

IDE 3 531.712 3.399 

NGR 3 531.710 3.397 



TABLE II: The comparisons among various cosmological models through the same method and observations. 



and more independent cosmic methods can give a more confirmed discrimination. We further compare NGR model 
with some of current successful dark energy models including Chevallier-Polarski-Linder (GPL) parametrization |34| . 



generalized Ghaplygin gas (GGG) and interacting dark energy (IDE) model [35!] through AIG information criterion. 
The AIG is defined as AIC = Xmin + ^/c, where k is the number of parameters. We show the comparisons in Table. |lll 
NGR model as a three-parameter cosmological model can compete with GPL and IDE model. From the discussions 
above, we can see NGR model gives a good discrimination between RDE model and XGDM model. Besides, as a 
unification of RDE model and XGDM model, it can be a good choice for describing the universe itself. 
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FIG. 6: The 2-D regions and 1-D marginalized distribution with the 1-a and 2-cr contours of parameters fimo and a in NGR 
model where we fix /3 = 0.5, for the data sets SNe+CMB+BAO. 
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